Frustration effect on magnetic excitations in a two-leg spin-ladder system 
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We theoretically study the magnetic excitations in a frustrated two-leg spin-ladder system, 
in which antiferromagnetic exchange interactions act on the nearest-neighbor and next-nearest- 
neighbor bonds in the leg direction, and on the nearest-neighbor bonds in the rung direction. 
A dynamical spin correlation function at zero temperature is calculated by using the dynamical 
density-matrix renormalization-group method for possible magnetic phases, i.e., columnar-dimcr 
and rung-singlet phases. The columnar-dimer phase is characterized by multi-spinon excitations 
with spin gap, while the rung-singlet phase is dominated by the triplet excitation in the rung direc- 
tion. It is found that a major difference between these two phases appears in the spectral weight of 
magnetic excitations, in particular, of the bonding and anti-bonding modes in the rung direction. 
Therefore, we can distinguish one phase from the other by the difference of the spectral weight. 
Furthermore, we examine frustration effect on both modes in the rung-singlet phase with a per- 
turbation theory from the strong coupling limit. It is shown that the anti-bonding mode is stable 
against the frustration, and a wavenumber with minimum excitation energy is shifted from the 
commensurate to incommensurate ones. On the other hand, the bonding mode is merged into the 
continuum excitation of multiple triplet excitations by increasing the frustration. By comparing our 
results with inelastic neutron scattering experiments for BiCu2PC>6, the magnitude of the magnetic 
exchange interactions and the ground state will be determined. 



I. INTRODUCTION 

Spinon is an elementary excitation characteristic to 
the strongly correlated electron systems in low dimen- 
sions-^, and was observed by angle-resolved photoemis- 
sion spectroscopy 2 - and neutron scattering 3 -. It is a man- 
ifestation of electron fractionalization called spin-charge 
separation, in which the electronic excitation is split into 
the spin degree of freedom (spinon) and the charge one 
(holon), instead of a quasiparticle. In the frustrated spin 
systems with exotic ground states such as spin-liquid and 
valence-bond-solid states, the spinon is associated with 
a topological excitation similar to the magnetic domain 
wall 4 . In contrast to the collective excitation in a mag- 
netic long-ranged ordered state, i.e., magnon, the spinon 
carries a fractional quantum number, spin 1/2, accompa- 
nied by a spin string. Because of this topological nature, 
two spinons make a bound state called triplon in some 
cases^. 

The frustrated two- leg spin-ladder system, i.e., J1-J2- 
Jp model (Fig. [T|), provides those exotic magnetic states, 
since this model bridges two different systems: the frus- 
trated spin chain, i.e., J1-J2 model and the unfrustratcd 
spin ladder, i.e., J\-J p model. Here, J\ and J 2 are the 
nearest and the next-nearest antiferromagnetic exchange 
interactions, respectively, and J p is the antiferromagnetic 
interaction in the rung direction. The ground-state phase 
diagram has been examined numerically^, and has two 
phases; the columnar-dimer (CD) phase, which continu- 
ously merges to the dimer phase of the frustrated chain 
system (J1-J2 model) in the zero limit of J p , and the 
rung-singlet (RS) phase, which continuously changes to 



the ground state of the unfrustrated ladder system ( J\- 
J p model) in the zero limit of J 2 - Note that the phase 
transition between the RS and the CD phases is of the 
Ising-type^. On the other hand, the excited states of the 
J\-J%-Jp model will be non-trivial, since the frustration 
leads to the confinement of spinons and triplons 4 -. The 
two limiting cases of the J1-J2 and J\-J p models are in- 
volved in the J\-J%- J p model. Therefore, this frustrated 
two-leg spin-ladder is useful for understanding the frus- 
tration effect on magnetic excitations. Several properties 
of the two limiting cases are summarized below. 

The J1-J2 model has two ground states: gapless spin- 
liquid and gapped dimer phases. Between these phases, 
Kosterlitz-Thouless transition occurs by changing the ra- 
tio of the exchange interactions, J 2 /Ji—~—. In the spin- 
liquid phase, spin correlation function exhibits a critical 
behavior and has no long-range order. The excitations 
are understood well with the spinon picture, which is ob- 
tained by the Jordan- Wigner transformation of the spin 
operators. In particular, the triplet excitation is domi- 
nated by two-spinon excitatio n 10 ^ 1 . In the dimer phase, 
neighboring two spins constitute a singlet dimer, whose 
dimer correlation remains finite in the limit of infinite 
distanc o 12 ' 13 . The ground state is doubly degenerated 
because of two choices to form a singlet dimer with ei- 
ther a right neighboring site or a left one. 

The intensive studies on the unfrustrated ladder sys- 
tem, i.e., the J\-J p model, show that one-dimensional 
spin-liquid phase is fragile with respect to the inter- 
chain interaction^. The RS phase appears in the ground 
stated, where two spins connecting a rung bond consti- 
tute a singlet dimer and there is no degeneracy. Thus, 
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the RS phase is another gapped phase different from the 
dimer phase. The magnetic excitation is dominated by 
triplon, which is a bosonic quasiparticle defined as the 
triplet excitations from the RS dimer a 15 ' 16 . It is also in- 
terpreted as a bound state of two spinons 3 - - — . The triplon 
always appears as an anti-bonding mode in the rung di- 
rection due to the parity changes between the ground 
state and the triplon excited stated. On the other hand, 
the excitation with even number of triplons appears as a 
bonding mode in the rung direction. The bound states 
of two triplons are important for low-energy excitation 
in the bonding modes, since the intensity is much larger 
than the two-triplon continuum^. 

In this paper, we investigate the magnetic excita- 
tions in the J\-J%-J v model. Our main purpose is to 
make clear the frustration effects by Ji on the excita- 
tion spectra, which are observed by inelastic neutron 
scattering for BiCu2P06 (Rcf. 18-21). The compound 
should be compared with the unfrustrated two-leg spin- 
ladder compound, SrCu2 03 (Ref. 22). Not only the 
dispersion relation but also the spectral weight is cru- 
cial for analyzing experimental data. Thus, we calcu- 
late dynamical spin correlation function (DSCF) by dy- 
namical density-matrix renormalization-group (dynami- 
cal DMRG) method. Furthermore, we propose how to 
determine the ground-state phase and how to estimate 
the interaction J\, J2 and J p in BiC^POg, by combining 
the DSCF and inelastic-neutron-scattering data. Our an- 
other purpose is to clarify difference between two phases, 
i.e., CD and RS phases, by the DSCF and to obtain re- 
lation of the exchange energies and the spectrum. 

The rest of the paper is organized as follows. First, 
we present the model Hamiltonian of the frustrated two- 
leg spin-ladder system and method to calculate DSCF 
by dynamical DMRG in Sec. [TTJ We also prepare some 
notations for a perturbation theory to explain behaviors 
of the excitations. In Sec. Mil we show the numerical re- 
sults and some interpretations on magnetic excitations. 
We discuss the frustration effects on the excitation. The 
ways to estimate the magnitude of the exchange inter- 
actions are also discussed. Finally, we summarize the 
results and discuss consistency with previous studies in 
Sec. El 



II. MODEL AND METHOD 

The model Hamiltonian of the frustrated two-leg spin- 
ladder system in Fig. Q]is defined as 
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FIG. 1. A frustrated two-leg spin-ladder system. Balls repre- 
sent spin-i sites. The nearest- and next-nearest-neighbor in- 
teractions in the leg direction, Ji and J2, are associated with 
solid and dashed lines, respectively. Double solid line denotes 
the nearest-neighbor interaction in the rung direction J p . 



where Ji(> 0) and h(> 0) are the magnitude of antifcr- 
romagnetic nearest-neighbor and next-nearest-neighbor 
exchange interactions, respectively, and J p (> 0) is that 
of antiferromagnetic nearest-neighbor interaction in the 
rung direction. Sj tU ^ is the S = 1/2 spin operator on j 
site in the upper (lower) chain. 

To investigate magnetic excitations in the frustrated 
ladder system, we calculate the z component of DSCF at 
zero temperature given by 

X (q,u>) = —3 / dte^(0\S z \q,t)S z (q,0)\0) 



= - lim *9<0|S«(-«l) ] . 

7-^0 it u — H + e + vy 



S z (q)\0), 
(5) 

where |0) is the ground state. The x (y) component 
equals to the z component because of no anisotropy in the 
model Hamiltonian ([T]). To obtain the excitation spectra, 
we set the infinitesimal value 7 to be 0.1, which is enough 
to examine the qualitative behaviors of the spin excita- 
tions in this system. The momentum representation of 
the spin operator with the open boundary condition is 
given by, 



S z (q) 
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with 



S lq v =0 ~ n( S j,u + S j,l)' 



S 3,qy=n - O^j.U S j,l)> 



(6) 

(7) 
(8) 



where x axis is set to the chain direction and y axis is in 
the rung direction. N is the rung number of the system, 
namely the total number of spin sites is 2N . 

We use the dynamical DMRG metho d 23 ' 24 to obtain 
the DSCF. This method requires three target states: |0), 
5 z (q)|0), and [w--H + e o + «7]" 1 5" z (g)|0). The last target 
state, so-called the correction vector, is obtained by using 
a modified conjugate gradient method in the dynamical 
DMRG algorithm. 
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We check the convergence of the ground-state energy 
for the DMRG truncation number m. We confirm that 
the ground-state energy per site converges within a nu- 
merical error less than 10~ 3 for m = 150 and use this 
value to obtain magnetic excitation spectra in a 32-rung 
frustrated ladder system. 

A perturbation analysis with resonated valence bond 
(RVB) picture is also used to obtain the ground state and 
low-energy excited states. We can rewrite the perturbed 
Hamiltonian, hj (R = 1,2), by the projection opera- 
tors onto RVB configurations. The RVB picture helps 
us to handle complicated spin-configurations easier (see 
Appendix). We prepare the notations for singlet pair, 
[j, k] p , and triplet one with S z = 0, {j, k} p , of two spins 
on j and k sites (j < k) as follows, 

b',fc] P =i,u,l = ;^=(| t>j,«J l)k,b p - | l)j,a p \ Dk,bp), (9) 
{j,fc}p=i,u,l = i(| t)j,a,| l)k,b p + I l)j,a p \ t)fe,6 P )(10) 

where the index p denotes the position of singlet or triplet 
bond, i.e., p — i indicates an inter-chain bond, and p = 
u(l) means a bond in the upper (lower) chain. The local 
spin state on j site in the upper (lower) chain is given by 
I t)j,u(l) or I i)j,u(i)- Then, a, = u and h = 1 for p = i, 
while a u = b u = u (oi = b\ = 1) for p = u(l). 

Below, the index p is omitted for the intra-chain, if it 
is obvious. We define two types of singlet configurations 
between j and k rungs as, 



(j,fc)r = \j,j]i[k,k]i, 
(j,k)\ = [j,k] u [j 7 k] h 



(11) 
(12) 



By using these expression, S z = states in total 5 = 
and 1 subspaces with two triplets on j and k rungs are 
obtained as, 



2 / 1 



(13) 



III. RESULTS 

In this section, we present numerical results of mag- 
netic excitations in the frustrated two-leg spin-ladder 
system, and discuss the frustration effects by Ji. Fig- 
ure [2] shows the DSCF in a 32-rung frustrated lad- 
der for three phases: the incommensurate CD phase 
(Jp/Ji = 0.2, J2/J1 = 0.6), the incommensurate RS 
phase {J p /J\ = 1.0, J2/J1 = 0.6), and the commensu- 
rate RS phase (J p /Ji = 1.0, J 2 /Ji = 0.1). 

First, we can find an obvious difference between the CD 
phase and the RS phase in the DSCF. That is, the spec- 
trum of the bonding mode in the CD phase (Fig. [2a)) 
exhibits almost the same shape with the same intensity 
as that of the anti-bonding mode (Fig.[2jb)), while those 
are completely different in the RS phase (Fig. [2Jc)-(f)). 
The difference between the CD phase and the RS phase 
is understood as the difference of the elementary exci- 
tations in both phases. In the CD phase, the key ele- 
mentary excitation is spinon, which strongly reflects one 
dimensionality. Thus, there is few difference between the 
bonding and anti-bonding modes in the CD phase. On 
the other hand, the major elementary excitation in the 
RS phase is triplon, which is the excitations from singlet 
to triplet in rung bonds. The parity of the triplet state is 
different from that of the singlet state. Therefore, the ex- 
citation exchanging the rung parity, that is anti-bonding 
mode (q y = 7r), is dominant in the RS phase. With fi- 
nite Ji, the ground state is mixed with singlet states in 
subspaces of even-number triplets, which has the same 
parity as the even-number triplet excitations. Hence, the 
excitation of the bonding mode (q y = 0) appears with 
finite J\. We stress that the difference in the spectral 
weights plays a key role in identifying the ground-state 
phase by using inclastic-ncutron-scattering data for the 
frustrated two-leg spin-ladder compound, BiC^POg. 

In the following subsections, we clarify the origin of 
the spectra for each phase with three excitations: spinon, 
triplon, and bound state of two triplons. The frustration 
effects are also clarified. 



and 



0*. fc )i = y| [V, k}u[j, % + [j, k]u{j, k}i j , (14) 

respectively. To simplify expressions of the perturbation 
theory, we use the Fourier transformation of one- and 
two-triplet excited state given by, 

|1 : q x ) ^E^fw'l'IlKH ( 15 ) 



and 



\2:S,{( k ,L)) = -L y £e™(j,j + L) s ]J [m,m]i 



(16) 

respectively, where S = 0, 1 and L = 1, 2, • ■ ■ , (N — l)/2. 



A. Multi-spinon excitations 

The phase transition between the CD and RS phases in 
the frustrated ladder system belongs to the same univer- 
sality class as the Ising transition^. The crucial point of 
the transition is a change of the ground-state degeneracy 
and the symmetry. Although the ground state has no de- 
generacy in the RS phase, there are twofold-degenerate 
ground states in the CD phase. At zero temperature, the 
translation symmetry is spontaneously broken in the CD 
phase. The CD phase continuously changes to the dimer 
phase in the J1-J2 model in the zero-coupling limit in the 
rung direction. Therefore, the key excitation should be 
spinon and is almost identical to the elementary excita- 
tion in the frustrated chain system. 

In Figs. [5] (a) and (b), we can see the lower edge of 
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FIG. 2. (Color online) The dynamical spin correlation functions (DSCF) of a 32-rung frustrated ladder for three phases: (a), (b) 
the incommensurate CD phase (J p /Ji = 0.2, J2/J1 = 0.6), (c), (d) the incommensurate RS phase (J p /Ji = 1.0, J2/J1 = 0.6), 
and (e), (f) the commensurate RS phase (J p /Ji = 1.0, J2/J1 = 0.1). (a), (c), and (e) [(b), (d), and (f)] represent the excitations 
with q y — [q v = tv^. We set the DMRG truncation number m = 150 and the broadening factor 7 = 0.1. 



the two-spinon continuum, which corresponds to the des 
Cloizcaux- Pearson mode in the unfrustrated limit J2 — >• 
in the chain syste m 11 ' 26 . The excitation has a minimum 
energy at an incommensurate wavenumber q* near the 
classical solution q*r ci \ = cos _1 (— Ji/4J 2 ) = q*. The 
dispersion relation is almost symmetric with respect to 
q x = 7r/2, although the intensity is weaker for q x < n/2 
than q x > tt/2. This feature originates from the dou- 
bled unit cell due to the dimerization. To explain this 
briefly, we consider the one-dimensional version of the 
DSCF for the frustrated chain system at the Majumdar- 
Ghosh point J2/J1 = 0.5 with the periodic boundary 
condition. The ground states doubly degenerate as fol- 
lows: 

N/2-1 

|MG),= [ [2.7,2.7 + 1], (17) 



3=0 
N/2-1 

\MG) r = [ [2j + l,2j 



2]- 



(18) 



Thus, we obtain the relation between x(li w)( 1D ) and 
x(tt - q,cj)( 1B 1 as, 



x(tt - q,u) 



(id) = 1 + c os(g) 
1 — cos(g) 



(ID) 



(21) 



where q =/= 0. For q = [7r/2,7r], we can easily confirm the 
relation about the spectral weight as, xi 71 ~ <z,w)' 1D ' < 
x(q,w)( 1D \ with the same poles. This is due to the fact 
that the doubled periodicity of the dimer state, which 
occurs spontaneously, folds the Brillouin zone in half. 
In other words, the original dispersion relation with- 
out dimerization and its inverted one with respect to 
q = [tt/2,tt] arc overlapped each other in the Brillouin 
zone. Figures [5] (a) and (b) have such a situation. 



We assume that the system size TV is even. The ground- 
state energy of the Majumdar- Ghosh Hamiltonian is eas- 
ily obtained as cmg = — § NJ\. The one-dimensional ver- 
sion of the DSCF for the left MG state |MG) Z is given 
by, 



with, 



(ID) 



1 — cos(g) 
2^ 



1 



(n 



MG 



£mg) + «0 

(19) 



N/2-1 



k) = E e 2 ^{2j,2 J + l}n[2m,2m+l]. 



(20) 
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B. Triplon excitation 



We discuss the anti-bonding mode (q y = it) in the RS 
phase (Fig.[2](d) and (f)) in terms of triplon. The triplon 
excitation is defined as the rung triplet excitations in this 
study. In the perturbation theory from the strong cou- 
pling limit in the rung direction, the dispersion relation 
of the triplon within the second order is obtained as 
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1 + ai cos(q x ) + a.2 cos(2q x 



\)+0(a 3 



(22) 



r 



where a.\ = Ji/Jp and a?2 = Jij J p . The wavenmnber 
at the minimum energy of the triplon excitation q* cor- 
responds to the incommensurability determined by the 
classical solution <3£( c n. 



In Fig. [2] (d), the wavenumber of the lowest excita- 
tion q* = (2/3 ± l/33)7r is obtained for the 32-rung lad- 
der system. The wavenumber is slightly shifted from the 
classical one q x ^ — 0.6377T. It is known that the quan- 
tum correction to the wavenumber can be obtained by a 
mean-field analysis with rung-bond operators. However, 
the analysis requires proper renormalization to obtain a 
precise energy, which is, however, demanding. Thus, the 
numerical calculation is indispensable to obtain the pre- 
cise wavenumber with the minimum energy in the disper- 
sion relation of the triplon, especially, for large oti and 
a.2 region. 

I 



In Fig. [5](d), we find that the intensity is concentrated 
at the wavenumber q ~ 27r/3 with the minimum energy. 
In Fig. [2] (f), both the maximum intensity and the mini- 
mum energy shift from q ~ 27r/3 to q = tt simultaneously. 
Thus, we conclude that the frustration affects only the 
wavenumber and the triplon picture is still stable. 



C. Bound state of two triplons 

In this subsection, we discuss the bonding mode {q y = 
0) in the RS phase (Figs. [2] (c) and (e)). First, we show 
that the bonding mode disappears in the strong coupling 
limit J p / Ji — > oo in the rung direction without frustra- 
tion J 2 = 0. The ground state in this limit is given by 
the direct product of the rung singlets, |0)(°) = Iljlj'jli- 
In addition, the excited state at q y = in Eq. (0, 
S z (q x ,q y = 0)|0)(°), is explicitly rewritten as 



S*(q x ,0)\0) iO) 



1 



3 



=o)([7, 



i]i n [m,m] 



(23) 



because Sj q \j, j]i = 0. Therefore, the intensity of the 
bonding mode should vanish in the strong coupling limit 
in the rung direction. 

Next, we consider the perturbation effect of u\ = 
Ji/Jp <C 1 on S = 1 excitation in the bonding mode. 
The ground state within the first order of the perturba- 
tion with respect to a% is obtained as, 



V3 



> (1) = -^aiV7V|2: 0,(0,1)). 



(24) 



The excited state from this state exists: S z (q x , O^O^ 1 ' 7^ 
0. For example, the excited state at q x = tt is obtained 
as, 

^(7r,0)|0)W = -^|2:l,(7r,l)). (25) 

This is the origin of the bonding mode, so-called bound 
triplon. Since the bound triplon appears proportional to 
the first order of cei, the DSCF of the bound triplon starts 
from the second order. Figure [3] shows the a% dependence 
of the bound triplon at q x = tt obtained by the dynamical 
DMRG calculation for a 32-rung unfrustrated ladder with 
the truncation number m > 300. In the inset of Fig. [3] 
(b), we can see that the peak energy is proportional to 
aj. 

We can obtain the excitation energy of the S = 1 state 
of two triplons at q x = tt by the perturbation theory. In 



addition to the bound triplon whose energy is £bt /Jp = 
2 — i«i + |af + 0(a\ ) (Fig. [3] (b)), there is a continuum 
of two triplons starting from the energy, e; / J p = 2 + 
|af +0(a\ ). In fact, we can see the continuum above the 
bound triplon (the inset of Fig. [3] (a)). We note that the 
bound triplon is stabilized because of the dimerization 
energy in the leg, -\j x + \ J X = -\J\. 



The frustration effect is also investigated with the per- 
turbation theory. Figure U] shows the «2 dependence of 
DSCF at q x = tt with ot\ = 0.5. We can see that the en- 
ergy of the bound triplon shifts to the continuum of two 
triplons as oti increases. At the same time, the weight 
of the bound triplon moves into the continuum. Within 
the second order of ol\ and in the first order of «2, i-e., 
oli/olx-^X is assumed, c'bt/Jp = 



Ct2 



0(al,al,aia 2 ) at q x 



— 2 2^1 ~ r 
■ tt by the perturbation analysis. 
In the inset of Fig. [H we can make sure that the energy 
of the bound triplon follows the equation. We expect the 
collective mode to be weakened and to be buried in the 
continuum with large J%/ J\ in the strong coupling limit 
in the rung direction. Therefore, the bound triplon is 
unstable with large frustration, while the triplon is still 
stable with a slight shift of the wavenumber. 
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FIG. 3. (Color online) (a) ot\ dependence of the DSCF in 
a 32-rung unfrustrated ladder, namely J2 = 0, at (q x , q v ) — 
(327r/33, 0) in the commensurate RS phase. The inset of (a) 
shows the DSCF for ai — 1/2 (red solid line) and Lorentzian 
fitting (green dotted line), (b) The energy of the peak versus 
ai where the energy of the bound triplon with second order 
perturbation theory, e^, / J p — 2 — ai/2 + 9qi / 8. The inset of 
(b) shows the DSCF at the peak versus ai. Blue line shows a 
fitting function proportional to a\ with the coefficient, 0.265± 
0.002. We set the DMRG truncation number m = 300 and 
the broadening factor 7 = 0.05 to distinguish the high-energy 
continuum from the main intensity of single Lorentzian. 



IV. SUMMARY AND DISCUSSION 

The frustrated two-leg spin-ladder system was investi- 
gated as a quantum spin system bridging two different 
systems: the frustrated chain system and the unfrus- 
trated ladder system. We examined the dynamical spin 
correlation function at zero temperature to understand 
the elementary excitations and the frustration effects, by 
using dynamical density-matrix rcnormalization-group 
method. We showed the excitation spectra for three 
phases: the incommensurate columnar-dimer phase, the 
incommensurate rung-singlet phase and the commensu- 
rate rung-singlet phase. 

First, we made clear that the q y dependence of the 
spectrum is different between the columnar-dimer phase 
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FIG. 4. (Color online) a 2 dependence of the DSCF in a 32- 
rung frustrated ladder at (q x ,q y ) = (327r/33, 0) at «i = 0.5. 
The inset shows the energy of the peak, corresponding to 
the energy of the bound triplon, as a function of Q2. Here, 
egrp' 1 ' = £ B t(qi) + Q2 J p (green line). We set the DMRG 
truncation number m = 300 and the broadening factor 7 = 
0.05 



and the rung-singlet phase. This is because the elemen- 
tary excitations are different between the two phases. A 
key excitation in the columnar-dimer phase is spinon, 
which is the same as the frustrated chain system. In 
contrast, a major excitation is triplon in the rung-singlet 
phase. The difference was explained by using both spinon 
picture for the columnar-dimer phase and triplon picture 
for the rung-singlet phase. With the help of the pertur- 
bation theory from the strong coupling limit in the rung 
direction, we clarified the instability of the bound triplon 
with the frustration, although the triplon is still stable. 

This is different from a frustration effect on the triplon 
excitations in the ladder system^, where the frustration 
is introduced with a diagonal antiferromagnetic interac- 
tion ( J<i) in the nearest-ncighbor-bond square. This sys- 
tem corresponds to a frustrated chain with alternating 
the nearest-neighbor interaction^—. The RS phase ap- 
pears without Jrf, and the dimer phase appears at J<i = J\ 
with large J\ in the frustrated ladder system. Therefore, 
from a viewpoint on the ground-state phase transition, 
the in-chain frustration ( J2) effect is similar to the diag- 
onal frustration ( Jd) effect for the ladder system. How- 
ever, the frustration effects on the excitations are dif- 
ferent, because the triplon excitation is fragile with the 
diagonal frustration, but not with the in-chain frustra- 
tion as shown in the present study. Thus, the frustration 
effects on the excitations are important for clarifying the 
influence of the in-chain frustration in the ladder system. 

Experimentally, this model will be useful for un- 
derstanding the low-energy physics in a newly synthe- 
sized compound BiCu2P06. The magnetic structure of 
BiCu2P06 is characterized by \ spins on Cu sites ar- 
ranged in zigzag chains in the b axis^— . The zigzag 
chains are stacked up in the c axis. The exchange inter- 
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actions J\ and Ji are denned as the nearest-neighbor and 
next-nearest-neighbor interactions in the zigzag chain. In 
addition, there are comparable inter-chain interactions 
in BiCu2PC>6. Since the distance between chains is al- 
ternate, we define the rung exchange interactions J p and 
Jp as the nearest-neighbor interactions on the shorter 
and longer bonds, respectively. The compound was re- 
ported to be similar to a unfrustratcd ladder system at 
first. That is, Ji and J4 are negligible as compared 
with J\. However, the following numerical studies have 
shown that the susceptibility can be understood well as 
that in a frustrated two-leg spin-ladder system, where 
the frustration is introduced as a next-nearest-neighbor 
antiferromagnetic interaction. Namely, the frustration, 
J2, should not be ignored. Besides, the band-structure 
calculation presented by Tsirlin et al£^ has claimed that 
J p interaction is comparable to J\ and J2 , although J p is 
rather negligible. By comparing our results and inelas- 
tic neutron scattering experiments^, the magnitude of 
magnetic exchange interactions and the ground state of 
such a compound will be determined in the near future. 
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with these equations. 



h) >{j,j}i\j + R,j + R]> 
hfHjJHj + R,^ 
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rung triplet unstable. In particular, Eq. (|A.4[) determines 
the dispersion relation of the triplon in Eq. (I2"2l) . 
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Appendix: Perturbation theory with RVB picture 

We describe the perturbation-theory analysis in the 
rung-singlet phase for the two-leg spin-ladder system. 
Because of the symmetries of the Hamiltonian, i.e., trans- 
lation in the leg direction and reflection of two legs, RVB 
picture is useful for analyzing low-energy states with per- 
turbation theory. Some simple relations helps us to cal- 
culate high-order states and the energies with perturba- 
tion theory. In this paper, we concentrate on the ground 
state, one-triplet excitations, and two-triplet excitations. 

First, wc show some equations regarding to the rung- 
singlct configuration. In the rung-singlct phase, the 
ground state is described as a linear combination of some 
lcg-singlct-pair states. In practice, the perturbed Hamil- 
tonian acting on the rung-singlet configuration generates 
a leg-singlet-pair state: 

hf ) (j, j + R) T = \ (j, 3 + R)r - (j, j + R)i, (A.l) 

where R = 1,2. In addition, the perturbed Hamiltonian 
can generate a long-range lcg-singlct pair and crossing 
lcg-singlct pairs as follows, 

(A.2) 
(A.3) 

I 

For the one-triplet excited state, the perturbed Hamil- 
tonian can only move the triplet to another rung without 
generating leg-singlet pairs: 

(A.4) 
(A.5) 

I 

On the other hand, the perturbed Hamiltonian acting 
on a leg-triplet state generates a longer-range leg-triplet 
state as follows, 



I 

\\j + RJ + R]i(j, k )i - \\ k i j +R)h 

\{i,j + R)i{j, *>i - \{i,k)i{j,j + R)i. 



I 

= \ \j,Mj + R,j + R}u 

= \{j + RJ + R}i(j, k)i-^{k, k}i{j,j + R) h 
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^ R) {j,Mj + R,k)i = l\j + R,j + R]i(j,k)i, 



(A.6) 
(A.7) 



This means that the bound state, i.e., a nearest- neighbor 
leg-triplet state, is no longer stable with any perturba- 
tion. With only nearest-neighbor leg interaction hj, off- 
diagonal parts of the perturbed Hamiltonian for the mo- 
mentum distribution of leg-triplet state are canceled out 



I 

at q = 7r. This is the reason why the bound triplon is sta- 
ble around q = ir. The next-nearest-ncighbor leg inter ac- 

(2) 

tion hj mixes the nearest-neighbor leg-triplet state with 
longer-range leg-triplet states. Thus, the bound triplon 
is smeared out with large frustration. 
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